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Abstract: We extract diffractive parton distribution functions (DPDFs) and diffractive 
structure functions from the most recent HI and ZEUS diffractive DIS data obtained 
by various methods. We consider Pomeron as an object with parton distribution func- 
tion, evolving according to the next-to-leading order (NLO) DGLAP equations within the 
framework of the 'Fixed Flavour Number Scheme' (FFNS). Having performed a global fit 
analysis, we achieve a very good description of all available measurements by introducing 
a new set of quark distribution form for the Pomeron. We predict longitudinal and charm 
proton diffractive structure function as well. Our results are compared with other analysis 
from the literature. 
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1. Introduction 

Structure functions are key ingredient for deriving parton distribution functions (PDFs) 
in nucleons, which is important for any process which involves colliding hadrons. These 
PDFs allow us to predict cross sections at particle colliders. For the success of the physics 
program having precise knowledge of PDFs is required [||, |2|, ^ . Strategies to extract the 
most precise PDFs to be used at the LHC, and questions on how to use future LHC data 
to further constrain the PDFs are discussed in the PDF4LHC Q forum. 

The HI and ZEUS collaborations presented their results on inclusive and various ex- 
clusive reactions, which is being actively studied by theorists and give access to a broader 
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understanding of proton structure. Although data-taking there has been stopped, new 
results continue to appear. 

One of their most important experimental results, working at a center of mass energy 
of about 300 GeV is observation of a significant fraction, around 15%, of large rapidity gap 
events in deep inelastic scattering (DIS) |, |, |, |, |0|, |n|, 0, |l|, |l|, |l5|, |T|, 0, |18 



l9|, [20|, PI, p^, 23, 24 1 . In this event, which is called diffraction, an exchanged photon of 



virtuality dissociates through its interaction with the proton at squared four momentum 
transfer t to produce a hadronic system X with mass Mx- 

There are four different theoretical approaches to analyze diffractive data, which are 



described and compared in Ref |25]. These are: the Pomeron Structure Function (PSF) 
model formulated in the framework of Regge phenomenology, the Bartels-Ellis-Kowalski- 
Wusthoff (BEKW) two gluon exchange dipole model p^, p7|, Bialas-Peschanski (BP) model 
based on the BFKL Pomeron approach and the saturation model of Golec-Biernat and 
Wusthoff (GBW) [^9[. These four frameworks are based on completely different theoretical 
concepts. The best description of all available measurements can be achieved with either 
the PSF based model or the BEKW approach |25]. In addition, the PSF approach works 



good also when fitted to the ZEUS-Mx data set, which is not the case for other models 



1 25, 30 1 . In Regge language, diffraction and the rapidity gaps which persist at high energy 



are associated with Pomeron exchange; the structure of the Pomeron could then be clarified 



p. 

The study of DESY experimental data based on the different methods provides influen- 
tial information that shall use to achieve the best precision possible in extracting diffractive 
parton distribution functions (DPDFs). The purpose of the present study which is based 
on the PSF approach is to discuss extensively the treatment of data sets, test the compat- 
ibility of the data obtained with various experimental methods and bring deeper insights 
into Pomeron functional form. We perform a NLO QCD analysis of the most available 
HI and ZEUS observables in MS scheme in the framework of Regge phenomenology and 
extract DPDFs. 

The outline of the paper is to give an introduction to the theoretical framework adopted 
for the diffractive events in Section |2[ Section ^ provides a new technique of parameteriza- 
tion and present light and heavy formulation of the structure function. Section Q describes 
methods of selecting diffraction at HERA and tackles the technical issue of compatibil- 
ity between different data sets. We discuss our global fitting procedure and the method 
of minimization has been applied along with in Section ||. Our results are presented in 
Section |6[ Section |7| contains conclusion and introduces a FORTRAN-code which is available. 

2. Theoretical framework of the diffractive events 
2.1 Diffractive cross sections 



The data are often presented in the form of a t-integrated reduced diffractive neutral current 
, defined via 



cross section (To?*-^\ defined via 



dxjp d/3 dQ2 xQ^ 



l-y+y->ja^('\xjp,f^,Q^) , (2.1) 
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or in terms of a diffractive structure function F^^^\x]p, P,Q'^) 



^^^(3) = pDis) _ ^pm) ^ (2.2) 

where 1+ = 1 + (1 — y)^. Due to the presence of the /Yj^ factor, the second term in ( |2.2| ) 
can be neglected anywhere but at very large y |3^, Since the kinematical variables are 
bound by sxy = Q^, the longitudinal structure function plays an important role in low-x 
physics. The effects of F^^^^ and F^^^^ are considered through their NLO dependence on 
DPDFs which will be shown in next section. 

2.2 QCD hard scattering factorization 



The proof that QCD hard scattering factorization can be applied to diffractive DIS |35| 
indicates that the cross section for the diffractive process can be considered in terms of 
convolutions of universal partonic cross sections a^^{x,Q'^), which are the same as those 
for DIS, with DPDFs /f as 

da^P^^''^ix,Q\x]p,t) = J2 fF{x,Q\xp,t) ^ da'\x,Q^) . (2.3) 

i 

Moreover, it is assumed that the shape of the DPDFs is not dependent to xjp and t and 



their normalisation is governed by Regge asymptotics [16, |3^. The mentioned assumption 
is required and compatible with the data fitted. The DPDFs can then be presented as a 
sum of two factorised contributions corresponding to Pomeron and sub-leading Reggeon: 

fPix, Q\xjp, t) = f,p,p{xjp, t) ■ /,(/3, Q2) + nn • f,a/p{^jp, t) ■ ff{(3, Q^) . (2.4) 

fiiP: Q"^) and ff^iP, Q^) are the partonic structure of Pomeron and Reggeon, respectively. 
Reggeon contribution contributes significantly only at low (3 and large xjp. To have detailed 
information on the flux factor refer to Refs. |16, |l^, |37|| . 



3. Extraction of parton densities in the Pomeron 

In comparison to the inclusive parton distribution, the DPDFs contain two additional 
variables (xjp, t). Since they do not influence evolution, we focus our attention on the term 
depending on /3 and Q^. 

3.1 Parametrization 

Due to the availability of new high precision measurements, phenomenological groups try 
to provide for precise understanding of the nucleon structure and its partonic content. 
Following our previous studies on extraction of precise polarized and unpolarized PDFs 
1 38, 43], we have motivation to tackle with the parameterization of DPDFs. 



A simple instruction is adopted in which the parton distributions of both the Pomeron 
and the Reggeon are parametrised in terms of non-perturbative input distributions as a 
function of z at Qq. 
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The structure of sub-leading trajectory ff^ are obtained from a parameterisation extracted 
from fits to pion structure function data with a free normahzation to be determined by the 
diffractive data. We test the dependence on the pion structure function by performing our 
fitting procedure base on the the GRV parametrisation instead on the Owens one I^Sj . 
No considerable deviation exists between both fits (similar result is reported in [p5|]). 

We parametrize the Pomeron, a light flavour singlet distribution ^(z), consisting of u, 
d and s quarks and anti-quarks with u = d = s = u = d = s, and a gluon distribution g{z) 
at an initial scale Qq=3 GeV^, such that 

zJ:{z,QI)=Aj:z^^ {l-zf^{l + D^z + E^z^^) , (3.1) 

zgiz,Ql) = Agz''^il-zf^e-T^, (3.2) 
with a = 0.01. Here z is the fractional momentum of Pomeron carried by the struck parton. 



The DPDFs can be extracted from fits to the data applying the DGLAP [46| evolution 
equations. We present a new set of quark singlet distribution form which will be discussed 
briefly latter. 

The gluon density is poorly constrained by the data, therefore it found to be insensitive 
to the Bg parameter. Performing a fit without the parameter Cg confirms the lack of 
sensitivity of gluon distribution to large z. Finally, we consider the parameterization of 
the gluon density identical to that used in the second analysis of [16|. Our motivation to 



apply this simple form of gluon distribution is that predictions based on it give the better 
description of diffractive dijet production cross sections measured in DIS in comparison 
with the form including Cg parameter 

We studied the value of a and found that it is not fundamental, e.g. a = 0.001 also 
results in an acceptable fit. Moreover, we performed the analysis with a = considering 
the behavior at high z by substituting (1 — z)^'^ . We did not get any satisfactory results. It 
means that although the exponential term is large in a region where no data are considered 
in the fit, the form of the parameterization for the gluon distribution at large z is important 



3.2 Light flavour contribution 

The Eqs. ( |2.1| ) and ( |2.2D are written in analogy with the way dPa'^^^'^'^ /dxBj dQ^ is related 
to the structure functions F2 and for inclusive DIS. Similarly, diffractive structure 
function can be described by 

F,(/3,Q2) = ^C'/(/3,^)® /,(/?, /i^) , (3.3) 

to all order in perturbation theory. Here denotes the factorization scale. As outlined 
earlier fj is considered in terms of the light fiavour singlet distribution Ti{z) and gluon 
distribution g{z). Generally speaking, structure function can be written like 

i^i(/3,g') = i^'*^'*(/3,Q') + i^'^'^'^^(/3,g^mi) . (3.4) 
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The flavour singlet contribution up to NLO is given by ||4^, ^ 



where i = 2,L and the Cj^g and Cj^g are the common NLO coefficient functions |49, 50|. 
3.3 Heavy flavour contribution 

The treatment of heavy quarks is something that nearly every group does slightly and it 
can lead to surprisingly different results for PDFs extracted . We consider the effects of a 
heavy quark within the framework of the so called FFNS where, beside the gluon, only the 
light quarks q = u,d,s are considered as 'intrinsic' genuine partons, i.e. massless partons 
within the nucleon, and heavy quarks h = c,b,t should not be included in the parton 
structure of the nucleon. They are always produced in the final state from the initial light 



quarks and gluons in near threshold region, i.e. Q ~ mj^ |5l|, 52|. However, even for very 



large values of Q^, ^ m?, these FFNS predictions are in remarkable agreement with 



DIS data g, 

The heavy structure functions are given through [55| 



Ft{P,Q') = Y,cff'-^{QVml)^f;iiQ') . (3.6) 
k 

Here, all quark flavours below nih are treated as zero-mass and one sums over k = 
g,u,u,d,d, ... up to nj flavours of light (massless quarks). The mass of heavy quark, 
nih, appears explicitly in the Wilson coefficients ^ ' , as indicated in Eq. (|3.6D . 

No simple analytic expressions can be given for the coefficient functions. To consider 



heavy contribution of F2^l we follow the standard manner applied in [56| 



Ffc(/3,Q^m2^ 



+ 



/3 ^ 



^\^2 f ,^2^ (0) 



Z L 



+ 



i=q,q 



+ei,/,(^,^^)(4i+C^ln 



(3.7) 



where k = 2,L and the upper boundary on the integration is given by Zmax = Q'^/iQ'^ + 
Am?). Further fi{x,iJ,'^), {i = g,q,q) denote the parton densities in the proton, ej^ and e| 
represent the charge squared of the heavy and light quarks respectively. The coefficient 
functions, represented by cj^\{V:0 ^^'^ '^^kii'^^O ^'^kii'H^O are fully described in 

|57, 58| up to NLO. Some progress at next-to-next-to-leading order (NNLO) can be found 



m 
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4. Data Analysis 



4.1 Methods of selecting diffraction at HERA 

Cross sections in deep inelastic difFractive scattering are not uniquely defined. Different 
methods exist to select diffractive events. These methods select samples which contain 
different fractions of proton dissociative events. Cross sections are not always given with 
corrections for proton dissociation. Three distinct methods of 

• Forward (Leading) Proton Spectrometer Method, 

• Large Rapidity Gap Method, 

• Mx Method, 

have been employed by the HI and ZEUS experiments, which select inclusive diffractive 
events of type ep — >■ eXY. The advantages and disadvantages of each method are explained 
in detailed in |^, ||]. 

A reasonable compatibility between these techniques and between HI and ZEUS results 
have been observed, which shows that there is no strong bias between these experimental 
techniques. As can be seen in several experimental references 32, 60, 63], there 



are tolerable agreement and common fundamental features between different data sets for 
much of the kinematic range. However, there are clear regions of disagreement especially 
between the Mx and LRG Methods H, ||. 

The full HERA data sample analysis is a powerful technique to achieve the best precision 
possible in extracting DPDFs. First steps are taken towards the combination of the HI and 
ZEUS results However, already at the present level, much can be done with existing 
data for the understanding of diffraction at HERA. Our strategy to study the HERA data 
will be discussed in Section | 

4.2 Data sets 

In this Section, we describe the different available HI and ZEUS data sets. Since the 
various data sets correspond to different ranges in the outgoing proton system mass. My, 
i.e. My = TUp in case of EPS (LPS) results and My < 2.3 GeV in case of the Mx 
results, additional factors are required before comparisons can be made. For all data and fit 



comparisons, all data are transported to the Hl-LRG-06 [16| measurement range My < 1.6 
GeV. The scaling factors, which has been obtained by corresponding experimental groups, 
are independent of /?, and xp [31, 32 1. As it was discussed in the experimental analysis. 



this is simply impossible to access any sensitivity of proton dissociation with respect to 
kinematics. To summarize, we give the HI and ZEUS available different data sets together 
with their scaling factors in the following. 

• measured by HI collaboration using the Large Rapidity Gap Method labeled Hl- 
LRG-06 [^], this is the default data set which is not corrected further. The analysis 
is restricted to the region y < 0.9; 
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aj? measured by HI collaboration with the leading final state proton detected in 
Forward Proton Spectrometer labeled Hl-FPS-06 [17|, data multiplied by the global 
factor 1.23. The analysis is restricted to the region y < 0.5; 

aj? measured by HI collaboration with the leading final state proton detected in 
Forward Proton Spectrometer labeled Hl-FPS-10 [18|, data multiplied by the global 



factor 1.20. The analysis is restricted to the region y < 0.7; 

measured by ZEUS collaboration with the leading final state proton detected 
in Leading Proton Spectrometer labeled ZEUS-LPS-04 |jl^, data multiplied by the 
global factor 1/0.75=1.33; 

cjjP measured by ZEUS collaboration with the leading final state proton detected 
in Leading Proton Spectrometer labeled ZEUS-LPS-09 data multiplied by the 
global factor 1/0.75=1.33. The analysis is restricted to the region y < 0.5; 

measured by ZEUS collaboration using the Large Rapidity Gap Method labeled 
ZEUS-LRG-09 [^], data multiplied by the global factor 1.05 ^. The analysis is 
restricted to the region y < 0.6; 

measured by ZEUS collaboration using the Mx Method labeled ZEUS-Mx-05 
[^] , data multiplied by the global factor 0.88 ^; 

measured by ZEUS cohaboration using the Mx Method labeled ZEUS-Mx-08 
[p^ ], data multiplied by the global factor 0.88; 

F^if^ and a^^^^ recently measured by HI collaboration using the LRG Method la- 



beled Hl-LRG-11 p3[ |, cross sections data with proton beam energy of £'^=460, 575 
and 920 GeV multiplied by the global factors of 0.97, 0.99, 0.97, respectively. The 
analysis is restricted to the region y < 0.85; 

4.3 Our strategy on data combination 

There are several issues on data analysing that motivated us to try out a deeper look in 
comparison with what has been explained up to now. 

• The first challenge is that Mx, used by ZEUS to separate diffractive from non- 
difFractive events, and LRG methods use the same data and thus they are strongly 



correlated O, 62] 



^In Ref. |20j, the ZEUS cross section is quoted and given in the table at the proton mass. This imphes 
subtracting the proton dissociation background by applying a constant scaling factor 0.75. To correct them 
to the My range of Hl-LRG-06, My < 1.6 GeV, one has to multiply the cross section, before the subtraction 
of the proton dissociation background, by the factor 0.91. Moreover, the ZEUS data are higher than HI by 
13% on average. Therefore the global factor 0.91 x (1-0.13) x (1/0.75) ~1.05 must be applied. 

^In Ref. j2^, Mx results are normalized to ZEUS-LRG-09 results by the factor of 0.83. Consequently, to 
transport Mx data to the Hl-LRG-06 measurement the factor of 0. 83x0. 91x (1-0.13) x (1/0. 75)~0. 88 must 
be applied. 



-7- 



• Although the leading-proton data previously suffered from low statistics and hence 
were unlikely to have much influence on the fit results, the high statistics of the 
present data make them competitive in precision with the result of the LRG method 

• The LRG results from HI and ZEUS are compatible in most of the kinematic region 
covered by measurements [|6l|] . 

Due to the above items, data based on the Mx method has not been considered in our 
global analysis. To study the influence and compatibility of the leading-proton data, we 
perform two scenarios: 1) based on the both LRG and FPS/LPS data and 2.) only LRG 
data. Comparison between two scenarios shows consistency of leading-proton data and 
confirms the importance of including them to increase the fit precision. 
The total data sets that we use in the present global analyses together with their type and 
(3,X]p and range are listed in Table ||. The fitted normalization Mi of the data sets 
included in the global fit, for each data set i are also shown in this table. 



Lable 


Data set 


/?- range 


X ip-range 


Q^-range 


^Data 


M 


Hl-LRG-06 
Hl-FPS-06 
Hl-FPS-10 
ZEUS-LPS-04 
ZEUS-LPS-09 
ZEUS-LRG-09 
Hl-LRG-11 
Hl-LRG-11 
Hl-LRG-11 
Hl-LRG-11 


a^' 

pS(3) 
2,L 

'^r,(v^=225,252) 
C'r,(V5=319) 


0.004-0.8 
0.02-0.7 
0.006-0.562 
0.007-0.48 
0.013-0.609 
0.025-0.795 
0.089-0.699 
0.089-0.699 
0.089-0.699 
0.089-0.699 


0.001-0.03 
0.0011-0.08 
0.0025-0.075 
0.0005-0.06 
0.0009-0.09 
0.0005-0.014 
0.0005-0.003 
0.0005-0.003 
0.0005-0.003 
0.0005-0.003 


8.5-1600 
10.7-24 
8.8-200 
13.5-39 
14-40 
8.5-225 
11.5-44 
11.5-44 
11.5-44 
11.5-44 


190 |ie 

40 |17 
100 |1^ 
27 [19 
42 [20 
155 |2C 
20 ||23 
25 [23 
12 [23 
10 |23 


)| 

'1 

)| 


0.9958 
1.0023 
1.0002 
1.0005 
0.9845 
1.0094 
0.9739 
0.9739 
0.9608 
1 


Total 










621 





Table 1: Overview of the published data points for /3 0.8, > 2 GeV and ^ 3.5 GeV^ 
together with the fitted normalization shifts TVi. The normalisation uncertainty cancels in the ratio 



5. QCD fit 

Before getting into the detailed procedure of our QCD fit, let us mention a few points. 
Following the treatment of |l^], in order to avoid regions which are most likely to be 
influenced by higher twist contributions or other problems with the chosen theoretical 
framework, only data in the range /3 ^ 0.8, Mx > 2 GeV and ^ 8.5 GeV^ are considered 
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in the fit. We found that although individual fitting procedure of some of the data sets 
works reasonably well at lower Q^, the quality of the global fit with all available data sets 
drops for < 8.5 GeV^. The effect of are considered through its relation to the NLO 
parton densities, such that no explicit cut on y is required. 

The different fitting groups use various values of the charm and bottom quark masses, 
which may change PDFs up to a couple of percent O. MSTW |66| have recently presented 



a detailed study of mass dependence in the PDFs. Here, the heavy quark masses are set 
to world average values of nic = 1.41 GeV, mi, = 4.50 GeV and mt = 175 GeV [^]. The 
strong coupling is set via ai^^((5o) ~ 0.306 for three flavours, which corresponds to the 
world average ai^\M'^) = 0.1184 for five flavours [^8|. Our analysis is performed using 
the QCD-PEGASUS program [||]. We work at NLO in the QCD evolution and take the 
renormalization and factorization scales to be equal (/xr = fip). An N-space evolution 



of Eq. (3.2) would require an accurate internal re-parametrization which facilitates a 



computation of the moments on the Mellin inversion contour. 
5.1 Detailed procedure 

We step in to the process of this project by performing a QCD fit under the same conditions 
and conventions as in fit B of Ref. [^] . Having considered the same functional form and 
cut scenario for 190 data points, we obtain almost identical result in input scale of 2.5 
GeV2. 

In a second step, we combine all the data sets with their normalization factors as defined 



in Section HJ. It leads to a xVdof= 689.3/610 = 1.13. 

Finally, we find that the exponential term plays negligible role in quark singlet distribution. 
Moreover, considering the additional factor of (1 + Dj^z + E^.z^^) provides flexibility to 
obtain a good description of the data. Thus, we vary S distribution to the functional form 
of Eq. ( |3.1D to improve our fitting procedure. This specific choice of parameterization 
reduces x^/dof value from 1.13 to 0.99. 

As already discussed in |jl^, |7^, if the number of parameters in describing parton 
densities is small, the as well as the fit values would be sensitive to the choice of the 
input scale. However, it is expected that the fit result to be independent of Qq. Thus, 
another advantage of using the new form of DPDF, introduced in this paper, is that the 
fit result is independent of varying the input scale in the reasonable range (In the order of 
1 GeV2). 

aip(O) is considered as a free parameter in the fit. Since the experimental treatments 
and cuts are different for each data set, we assign a global Reggeon normalization parameter 
to each data set for which this contribution is required, similar method is applied in Ref. 
|25]. Considering Reggeon normalization parameter to Hl-LRG-11 is not needed because 



of its low xjp coverage. This gives us a total number of 14 unknown parameters which are 
presented in Table 0. 

Our result for aip(O) is compatible with the pomeron intercept describing soft hadronic 



scattering, ajp(O) ~ 1.08 |71, |7^, It is also consistent with that obtained from HI data 



previously measured using the LRG and FPS methods |16, 17, 18 1 and with the ZEUS 
measurements |21, 22 1. 
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Considering the general form of Eq. ( |3.lD for non-diffractive PDFs, in non-singlet 
combinations, e.g. the valence quarks and (u — d), B parameter is expected to be ~ 0.5. In 
singlet contributions, e.g. the sea and gluon, B is expected to be ~ [Q]. Here, in analogy 
with non-diffractive PDFs, we achieved similar amount for By,- 

The diffractive quark singlet and gluon distributions from our model are compared 
with the results from HI 2006 Fit B H, MRW |7| and ZEUS SJ H for different values 
of on a logarithmic z scale in Figure ||. At low Q^, both the quark singlet and the gluon 
densities remain large up to the highest z values accessed. 

As it was discussed, the contribution of the sub-leading Reggeon trajectories should 
be considered for xjp values substantially larger than 0.01. The contribution from these 
trajectories is modelled using the pion structure function. The pion PDFs are used in a 
region of low /? where they are not directly constrained by data. There is potential interfer- 
ence between the contributions from the Pomeron and sub-leading Regge trajectories with 
vacuum quantum numbers (e.g. the f-meson, see Eq. (14) of Ref. Q). In order to limit the 
influence of sub- leading Reggeon trajectory, we exclude data points with xjp > 0.01 and use 
a single A^jr for all required data sets. The results are shown in Figure ^. Only very small 
differences are observed between both fits. This analysis leads to ajp = 1.110 ± 0.0041. 



Parameters 


TKT 


As 


0.17 ±0.009 


By 


0.08 ±0.031 


Cs 


0.53 ±0.025 


Dy 


4.88 ±0.14 


Ey. 


-2.36 ±0.064 


Fy 


0.30 ±0.012 




0.44 ± 0.020 


aip{0) 


1.108 ±0.0035 


ArK(Hl-LRG-06) 


(1.22 ±0.18) X 10-'^ 


A^iR(Hl-FPS-06) 


(1.21 ±0.23) X 10~3 


A^iR(Hl-FPS-lO) 


(1.36 ±0.09) X 10-3 


AfiR(ZEUS-LPS-04) 


(1.64 ±0.26) X 10-3 


AfK(ZEUS-LPS-09) 


(2.25 ±0.12) X 10-3 


AfjR(ZEUS-LRG-09) 


(2.19 ±0.29) X 10^3 


xVdof 


601.92/607 = 0.99 



Table 2: Pomeron quark and gluon densities parameters and their statistical errors for combined 
data sets in MS scheme at the input scale Qq=3 GeV^. 
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5.2 The method of the minimization and error calculation 



The quality of fit is traditionally determined by the of the fit to the data |^5|, which is 
minimized using the MINUIT package . Xgiobai defined by 



Xgiobai(p) — E 



i=l 




T-\D,theor / \ ' 



D,data 



(5.1) 



where p denotes the set of independent parameters in the fit and n*^"*" is the number of 
data points included. For the i^^ experiment, p^j'^"-^"-^ ^-pD^ata pD,theor ^jgj^Q^g ^^le 
data value, measurement uncertainty and theoretical value for n*'^ data point. ATVj is the 
experimental normalization uncertainty and Afi is an overall normalization factor for the 
data of experiment i. We allow for a relative normalization shift Mi between different data 
sets within uncertainties AA/i quoted by the experiments. 

The errors include systematic and statistical uncertainties, being the total experimental 
error evaluated in quadrature. We check the fit stability by performing the two approaches 
with statistical and systematics errors added in quadrature or with statistical errors only. 
For a given set of data, the results based on the fit with statistical or total errors are 
very close. When moving to a combined fit of all data sets, although DPDFs show small 
differences between both fits, using statistical errors lead to fit with a large x^- 

There are clear procedures for propagating error experimental uncertainties on the 
fitted data points through to the PDF uncertainties. The most common is the Hessian 
approach. In this case we can consider 



global 



Xglobal 



Amin 



,{0) 



)(a,- 



(5.2) 



where the Hessian matrix is defined as 

Hij = 



1 



global 



2 daidaj 

The standard formula for linear error propagation is 

{AFf 



(5.3) 



(5.4) 



Since the derivative of F with respect to each parameter a, is required, this formula is 
not easily calculable. It can be improved by finding and rescaling the eigen vectors of H 
1 77, 78, ^]. In term of the rescaled eigevectors Zi, the increase in is given by 



2 _ 2 

Xglobal Xmin 



(5.5) 



The uncertainty on a quantity is then obtained applying 

iAFf = lY.[^(^t)-FiS-f^ 



(5.6) 
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Lable 




Data points 


Hl-LRG-06 

Hl-FPS-10 
ZEUS-LPS-04 
ZEUS-LPS-09 
ZEUS-LRG-09 

Hl-LRG-11 


0.79 

0.91 
0.23 
0.78 
1.08 
0.89 


190 iiel 

40 1171 
100 |18| 
27 |19l 
42 1^ 
155 |20| 
67 Q 



Table 3: values per data set for the gobal QCD fit witli statistical and systematic errors added 
in quadrature. 

where Sf and are PDF sets displaced along eigenvector directions by the given Ax^. 
The uncertainties on our DPDFs following this method are presented in Figure |l[ 

We present values for individual data sets in Table ^, which would allow some 
assessment of the degree of compatibility. It justifies our approach to combine all data 
sets. These fit results are displayed in Figure |3|. We show the quark and gluon densities in 
the Pomeron for individual HI and ZEUS data sets. We note that the Pomeron is gluon 
dominate for all fits. 

6. Discussion of fit results 

6.1 Behaviour of cross section and structure function 

The data on and F2 ■, which describes the structure of the Pomeron exchanged in the 
t-channel in diffraction, have two prominent characteristics p0[ : 

• Treatment versus /3 

As can be seen in Figures ^ and |5|, diffractive cross section is widely smooth in the 
measured /3 range. Considering the similarity between /3 in diffractive DIS and xsj 
in inclusive DIS, this is very different from the treatment of 0"^. Typical cross section 
effectively reduces for xsj ^ 0.2. 

• Treatment versus 

As illustrated in Figure |6|, the structure function F2 rises with for all /3 values 
(except the highest). This brings the scaling violations of F2 to mind, except that 
F2 increases with only for xsj ^ 0.2. In the proton, negative scaling violations 
manifest the existence of the valence quarks, though positive scaling violations are 
due to the growth of the sea quark and gluon densities. Consequently, the F2 data 
imply that the partons resolved in diffractive events are mostly gluons. 
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These results were already seen in HI and ZEUS papers in the mid 90's. Figures |7[ |8| and 
m display the reduced diffractive cross section as a function of for different regions of /3 
and xjp. Our model describes all the data well. 



Using the fit results of TKT superimposed on Mx data in Figures 10 and 11 shows 
that the fit can describe most kinematic regions. It leads to the conclusion that there seems 
to be compatibility between all data sets. 

6.2 Longitudinal structure function 

Predictions for are mainly determined by the form of the gluon density extracted from 
the fit. Consequently, the longitudinal structure function provides a way of studying the 
gluon distribution and a test of perturbative QCD This is considerably important 

at the lowest z values. 

Measurements of became possible following the reduced proton beam energy runs at 
the end of HERA operation 



and first results are recently presented in Ref. 



The reduced diffractive cross section as a function of /3 for the different proton beam 



energies [23] together with the prediction of our model are shown in Figure 12. Discrep' 



ancles of the measured cross sections from our F2 curves are considerable in the low and 
medium energy data. 



The measurements of F^ 



3[ are shown as a function of /3 in Figure 13. The data are 



compared with the predictions of our model and HI DPDF Fits A and B |16]. All three 
models are compatible with the data. 

The relative sizes of the diffractive cross sections, = F^/{F^ - ), is defined 



recently [23] by similarity to the inclusive DIS case 184]. The measurement of is shown 
as a function of (3 in Figure The data are consistent with our model. 

6.3 Charm contribution to the structure function 

Charm production in diffractive DIS has been measured by the HI and ZEUS collabora- 
tions. Due to the sensitivity of charm production to gluon- initiated processes this process 
is very important [^]. Charm quarks are selected by two independent methods in HI 
They are selected by the full reconstruction of D* mesons or by reconstructing the 



displacement of track, while just the later is used in ZEUS ]87]. In Figure 15 we present 
our results for charm diffractive cross section and structure function together with HI 2006 
DPDF l|l|] Fit A, Fit B and MRW |7|] and we compare them with HI and ZEUS data. 
ZEUS data and MRW curves are corrected with a factor of 1.23 to account for the differ- 
ence in the measured range from My = rap to My < 1.6 GeV ]17]. The predictions are in 
fair agreement with the data. 



7. Conclusion 

We have shown that the diffractive observables measured in the HI and ZEUS experiments 
at HERA can be well described by a perturbative QCD analysis which fundamental quark 
and gluon distributions, evolving according to the NLO DGLAP equations in FFNS, are 
assigned to the Pomeron and Reggeon exchanges. This work provides a detailed picture 
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of the Pomeron structure. In particular, a global analysis of all available data has been 
performed by introducing a new set of quark distribution form for the Pomeron. The new 
functional form developed here contains additional parameters provides flexibility to obtain 
a proper description of the data. We know that to constrain the gluon, the dijet data is 



required 88, ^ and this work is in progress. Having extracted the diffractive PDFs, we 
compute various diffractive structure functions. In general, we find good agreement with 
the experimental data, and our results are in accord with other determinations from the 
literature; collectively, this demonstrates progress of the field toward a detailed description 
of the Pomeron structure. 

Although these data obtained by various methods with very different systematics, they 
are broadly consistent in the shapes of the distribution throughout most of the phase space. 
This is a very important message from HERA that DPDFs are well compatible for both 
experiments. 

The new measurements performed by HI in a more extended kinematic regime will 
allow to further refine the results. Additionally, the Higgs boson may be produced at the 
LHC via a diffractive process in which fast protons are detected. A deeper understanding 
of diffraction, driven by the HERA result, could therefore aid in the discovery of the Higgs 
boson [||, 

A FORTRAN package containing our diffractive light and heavy structure functions 
F^^^^^^^^"'^^\ -^L*(iP*jRP^^^ ^^^^ Pomeron densities S and g with their errors at 

NLO in the MS scheme can be found in http://particles.ipm.ir/links/QCD.htni or 
obtained via e-mail from the authors. These functions are interpolated using cubic splines 
in and a linear interpolation in log (Q^). The package includes an example program to 
illustrate the use of the routines. 
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Figure 1: Comparison between ttie total quark singlet and gluon distributions obtained from our 
model (solid curve), HI 2006 DPDF Fit B (dashed curve) ZEUS SJ (dashesd-dotted curve) 
||l3|] and MRW (dashesd-dotted-dotted curve) |Q at four different values of as a function of z. 
The ZEUS SJ and MRW DPDFs plotted here are normalised to My < 1.6 GeV by muhiplying by 
a factor 1.23 relative to My = '7i„. 
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Figure 2: Comparison between the total quark singlet and gluon distributions obtained from our 
model. Results are presented with no cut on xp (solid) and xjp < 0.01 (dashed). 
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Figure 3: Singlet and gluon distributions of Pomeron as a function of z derived from QCD fits on 
Hl-FPS-10 m data alone, Hl-LRG-06 @ data alone, ZEUS-LRG-09 |20| data alone and all the 
data sets together. 
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Figure 4: The ZEUS-LPS-09 |^ diffractive cross section multiplied by xp, X]pa^^^\ as a function 
of /3 for different regions of and xp. The error bars indicate the statistical and systematic errors 
added in quadrature. The curves show our model reduced by a global factorl.33 to correct for the 



contributions of proton dissociation processes as described in Section 4.2 
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Figure 5: The Hl-LRG-06 reduced diffractive cross section, multiplied by a;jp, X]pa^^'^\ as 
a function of /? for different regions of and xp. The error bars indicate the statistical and 
systematic errors added in quadrature. The curves sho'w our model. The contribution of the 
sub-leading exchange alone is also sho'wn in xp = 0.03 and 0.01. 
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Figure 6: The ZEUS-LPS-04 |l9[ reduced diffractive structure function multiplied by xp, 
as a function of for different regions of /3 and xp. The error bars indicate the 
statistical and systematic errors added in quadrature. The curves show our model reduced by a 
global factor 1.33 to correct for the contributions of proton dissociation processes as described in 



Section 4.2 
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Figure 7: Comparison between the Hl-LRG-06 Q and ZEUS-LRG-09 Q measurements after 
correcting the latter data set to My < 1.6 GeV by applying scale factor of 1.05. The measurements 
are compared with our model (solid curve) and the results of the HI 2006 DPDF Fit B (dashed 
curve) [E6|, which was based on the HI data shown. 
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Figure 8: The Hl-FPS-10 ||l^ reduced difFractive cross section, multiplied by xp, X]pa^^'^\ as 
a lunction of for different regions of /? and xp. The error bars indicate the statistical and 
systematic errors added in quadrature. The curves show our model reduced by a global factor 1.20 



to correct for the contributions of proton dissociation processes as described in Section 4.2 
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Figure 9: The Hl-FPS-06 ||T^ reduced diffractive cross section, multiplied by xp, xpa^'^^\ as 
a function of for different regions of /3 and xp. The error bars indicate the statistical and 
systematic errors added in quadrature. The curves show our model reduced by a global factor 1.23 



to correct for the contributions of proton dissociation processes as described in Section 4.2 
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Figure 10: The ZEUS-A/j;-08 reduced diffractive structure function multiplied by xp, 
as a function of for different regions of /3 and xp. The error bars indicate the 
statistical and systematic errors added in quadrature. The curves show our model reduced by a 
global factor 0.88 to correct for the contributions of proton dissociation processes as described in 



Section 4.2 
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Figure 11: The ZEUS-Mx-05 diffractive structure function multiplied by xp, xjpFj as 
a function of /? for different regions of and xp. The error bars indicate the statistical and 
systematic errors added in quadrature. The curves show our model reduced by a global factor 0.88 



to correct for the contributions of proton dissociation processes as described in Section 4.2 
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Figure 12: The reduced diffractive cross section, multiplied by xp, xpa^'^^\ as a function of /? at 
fixed Q2 and xp for (from left to right) the 460, 575, 820 and 920 GeV data sets |2|. The error 
bars indicate the statistical and systematic errors added in quadrature. The curves show our model 
reduced by a global factor 0.97, 0.99 and 0.97 for Ep=A60, 575 and 920 GeV, respectively. 
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Figure 13: xpF^ measurement as a function of /3 measured at fixed and xp. The 
present fit is the solid curve. Also shown are the results of HI 2006 Fit A (dashed) and and Fit B 
(dashed-dotted) ||. 
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Figure 14: The ratio measurement i?^, as a function of (i at the indicated values of a; ip and 
. The error bars shown are the statistical and systematic uncertainties added in quadrature. Our 
result is the solid curve. 
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Figure 15: Comparison of our result (solid curve) for the contribution of the charm quarks to the 
diffractive cross section and structure function with HI 2006 DPDF n& Fit A (dashed curve), Fit 
B (dashed-dotted curve) and MRW (dashed-dotted-dotted curve) ]74| shown as a function of /3 for 
different values of Xip. The data obtained from the HI |86 and ZEUS D* production and from 
HI displaced track method [|6|. The error bars of the data points represent the statistical and 
systematic error in quadrature. Measurements at the same values of (3 are displaced for visibility. 
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